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ABSTRACT 


BACKGROUND AND OBJECTIVES: In this paper, a life insurance product is designed with 
the help of stochastic control approach. These products are defined in such a way that in 
exchange for receiving an amount as insurance premium that is paid at specified times, 
the insurer undertakes to pay insurance benefits when the insured is alive at the end of 
the contract. 

(METHODS: This research is an analytical study in terms of developmental-applicative 
purpose. In the literature of life insurance, there are various products that are not the 
same in the type of benefit payment and the timing of their implementation. Among 
these examples, term life insurance, term life insurance, and mixed life insurance can 
be mentioned. Traditional insurance products with fixed benefits are quickly losing 
their appeal due to inflationary markets. In this research, it is focused on the design of 
a life insurance product on the condition of life, which is connected to the investment 
markets. Stochastic differential calculus models have been used to simulate capital 
markets assets. All the numerical results of this research have been calculated with the 
help of Matlab and Maple software. 

FINDINGS: To achieve the best choice of investment type, with the help of stochastic 
optimal control tool, the best investment strategy was calculated for a person who has 
the CRRA utility function and buys this product, so that the most benefits are paid at the 
end of the contract. To invest in this contract, modeling was done in a non-risky market 
such as a bank and a risky market such as stocks, which have price jumps. In addition, 
to model the risk asset, the Merton model, which is a representative of the models with 
finite activity, was used, and at the end, a comparison was made for several mortality 
functions. 

CONCLUSION: The main purpose of this article is investment for the insured who bought 
this product. In the product designed in this article, the insurer undertakes to pay the 
premiums received at a guaranteed rate at the end of the contract. Also, the insured will 
share the profit from the investment based on a certain percentage that is determined 
at the beginning of each year. The simulations show that the behavior of the optimal 
consumption rate is the same as the Merton model with the approach that the behavior 
of full price jumps is transparent in the optimal consumption rate designed in this article. 
Investment results for several mortality functions are reported in the Numerical Results 
section. 


This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). 
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Fig. 1: The behavior of the optimal consumption rate for the De Moir mortality model with power utility function. 
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Fig. 2: The behavior of the optimal consumption rate for the De Moir mortality model with logarithm utility function. 


Optimal Consumption (c ) 


0.02 S| 


ice. Sie 

0 1 1 foe = ——_—— = — 

0 2 4 6 8 10 12 14 16 18 20 
Time(t) 


C=1.01, B=0.01 clr jell & 53 prolF prog Spo Joe Gly et Cengllrs eb b aingy Byers Ei LB) jlogei 
Fig. 3: The behavior of the optimal consumption rate with logarithm utility function for the Gampertz mortality model with 
B=0.01 and C=1.01. 
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Fig. 4: The behavior of the optimal consumption rate with power utility function for the Gampertz mortality model with B=0.01 
and C=1.01. 
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Table 1: The optimal investment strategy for the power utility function with different risk aversion coefficients. 
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g=3.2; csi=0.3; mu=0.28; sig=0.16; r=0.04; w=100; 
x=60; xx=20; pi=0.3901; T=20; fiVG=-0.1;; a=0.05; 
gg=2.75; fiKOU=0.7; 

t=1:1:T; 
si=(1-g)*(r+pi*(mu+t(0.5*(sig%*2))-r)- 
(0.5*g*((pi*sig)*2)))+fiKOU; 

for i=1:length(t) 

X=exp(csi*(w-x-T)-(si*(T-t(i)))/g); 

Y=(csi*(-g-1))* ((2*(w-x))/(w-x-t(i)))(1/g); 
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t(i)),g+1); 


XX=exp(csi*(w-xx-T)-(si*(T-t(i)))/gg); 
YY=(csi*(-gg-1))*((2*(w-xx))/(w-xx-t(i)))*(1/gg); 
GG=gammainc(csi*(w-xx),gg+1)-gammainc(csi*(w- 
xx-t(i)),gg+1); 


cc(i,:)=exp(a*(T-t(i))/gg)/(GG*XX*YY); 
c(i,:)=exp(a*(T-t(i))/g)/(G*X*Y); 

end 
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Consumption for Kou Jump 


g=8.04; csi=0.3; mu=0.28; sig=0.16; r=0.04; C=1.01; 
B=0.01; 

x=60; xx=30; pi=0.3901; T=20; fiVG=-0.01; a=0.05; 
gg=7.15; fiKOU=0.07; 

ZaribKou=[0.3463 0.7143 1.1053 1.5209 1.9626 
2.4321 2.9312 3.4617 4.0256 4.6252 5.2627 5.9405 
6.6614 7.4279 8.2431 9.1102 10.0325 11.0136 12.0574 
13.1678]; 

t=1:1:T; 
si=(1-g)*(r+pi*(mu+(0.5*(sig*%2))-r)- 
(0.5*g*((pi*sig)*2)))+fiKOU; 

for i=1:length(t) 


X=ZaribKou(i)* exp((1/g)*(((B/log(C))*(C*(x+t(i))- 
CAx))-si*(T-t(i)))); 
XX=ZaribKou(i)*exp((1/gg)*(((B/ 
log(C))*(CA(xx+t(i))-C4xx))-si* (T-t(i)))); 


ce(i,:)=exp(a*(T-t(i))/gg)/(XX); 
c(i,:)=exp(a*(T-t(i))/g)/(X); 
end 

hold on 

plot(t,c,.-’) 

plot(t,cc,’-+’) 
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end 

hold on 

plot(t,c,’-’) 

plot(t,cc,’-+’) 

ylabel(‘Optimal Consumption (c4*)’) 

xlabel(‘Time(t)’) 

AHWW“%“%“KLLVVVVVVLVLVVVLUHVVUV%W%%% 

AHWW““%““KKVVVVVVVVVVUVVVVWVVHV%VM% 

HWW%%% 

4 druloes (59) Cou! od Arbors digs CPL rel pj UF yo 

cul dg>g0 Maple jljiles 9 aS op F JS claysg) SoS 
lew] oS colin! 


> restart: 
s with( Student| NumericalAnalysis |) : with(LinearAlgebra) : 


o:= 0.16: 
p:= 0.5: 
A= 8.8: 
N:=0: 
uu := 0.28: 
r:= 0.04: 


gi=1:6'= 88:4 =2: 


>P0:= 0.5: 


fae (Spl te(P—1)) % 1 — 9842-2 —2-(1 —p) 1 — Pt) 8 — 9" 
-1) : 
3 
M2 = expand(Quadrature(f(t), t= 0 ..1, method = gaussian[3] , partition = 1, output 


=value)) : 
> 


M22 = expand(Quadrature(f(t), t= 0..1,method = newtoncotes|3]_, partition = 1, output 


=value)) : 
> 


> equ = P>—g-0"-P + M2: 
5 Optimal Strategy = (u + 0.5: os r) + Newton(equ(P), P= PO, tolerance = 10°) 


for i=1:length(t) 
X=mfun(‘Ei’,(B*(Cx))/log(C)); 
Y=mfun(‘Ei’,(B*(C4(x+t(i))))/log(C)); 
Z=exp(a*(T-t(i))-(B*C4(x+t(i)))/log(C))*log(C); 


XX=mfun(‘Ei’,(B*(C“xx))/log(C)); 
YY=mfun(‘Ei’,(B*(C4(xx+t(i))))/log(C)); 
ZZ=exp(a*(T-t(i))-(B*C%(xx+t(i)))/log(C))*log(C); 


c(i,:)=Z/(Y-X); 

cc(i,:)=ZZ/(YY-XX); 

end 

hold on 

plot(t,c,.-’) 

plot(t,cc,’-+’) 

ylabel(‘Optimal Consumption (c4*)’) 
xlabel(‘Time(t)’) 
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alpha=0.05; x=60; w=100; xx=30; 

t=1:1:20; 

for i=1:length(t) 
c(i,:)=((w-x)*(alpha*2))/((1+(alpha*t(i))-alpha*(w- 
x))+exp(alpha*t(i))*(1+alpha*(w-x))); 


cc(i,:)=((w-xx)*(alpha’2))/((1+(alpha*t(i))- 
alpha*(w-xx))+exp(alpha*t(i))*(1+alpha* (w-xx))); 
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